home *** CD-ROM | disk | FTP | other *** search
/ WINMX Assorted Textfiles / Ebooks.tar / Text - Mathematics - Numerical Mathematics and Computing (F).zip / xsinx.f < prev   
Text File  |  2002-06-11  |  773b  |  38 lines
  1. C
  2. C PAGE 63-65: NUMERICAL MATHEMATICS AND COMPUTING, CHENEY/KINCAID, 1985
  3. C
  4. C FILE: XSINX.FOR
  5. C
  6. C EXAMPLE OF PROGRAMMING F(X) = X - SIN(X) CAREFULLY (F)
  7. C
  8.       DOUBLE PRECISION DX,DY,DZ       
  9.       DX = 1.0D0/15.0D0     
  10.       DY = DSIN(DX) 
  11.       DZ = DX - DY
  12.       PRINT *,DX,DY,DZ      
  13.       X = 16.0    
  14.       DO 2 N = 1,52 
  15.         X = 0.25*X
  16.         Y = F(X)  
  17.         PRINT *,N,X,Y       
  18.    2  CONTINUE
  19.       X = SIN(1.0)
  20.       Y = F(1.0)
  21.       PRINT *,X,Y 
  22.       STOP
  23.       END 
  24.   
  25.       FUNCTION F(X) 
  26.       IF(ABS(X) .GE. 1.9) THEN
  27.         F = X - SIN(X)
  28.       ELSE
  29.         T = X**3/6.0
  30.         F = T     
  31.         DO 2 N = 1,9
  32.           T = -T*X*X/REAL((2*N+2)*(2*N+3))      
  33.           F = F + T 
  34.    2    CONTINUE  
  35.       END IF      
  36.       RETURN
  37.       END 
  38.